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Abstract 

Using the method of non-equihbrium statistical Zubarev operator an approach to the description 
of kinetics taking into account the nonhnear hydrodynamic fluctuations for quantum Bose system 
is proposed. The non-equilibrium statistical operator that consistently describes both the kinetic 
and nonlinear hydrodynamic fluctuations in quantum liquid is calculated. A kinetic equation for 
the non-equilibrium one-particle distribution function and generalized Fokker-Planck equation for 
non-equilibrium distribution function of hydrodynamic variables: densities of momentum, energy 
and particle numbers are obtained. A structural function of hydrodynamic fluctuations in cumulant 
representation is calculated. It provides a possibility to analyzing the generalized Fokker-Planck 
equation in Gaussian and higher approximations for dynamic correlations of hydrodynamic variables 
that is important to describe the quantum turbulent processes. The generalized Gross-Pitaevkii 
equation for the order parameter which takes into account nonlinear hydrodynamic fluctuations is 
obtained. 

PACS numbers: 67.40.-w, 47.37.+q 
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I. INTRODUCTION 



The examination of dynamic properties of quantum liquids as well as of the features 
of transition processes from gaseous state to fluid and superfluid one with temperature 
decrease remains a hard problem in modern physics. The development of the non-equilibrium 
statistical theory which would take into account one-particle and collective physical processes 
that occur in a system is an example of such problem. The separation of contributions 
from the kinetic and hydrodynamic fluctuations into time correlation functions, excitations 
spectrum, transport coefficients allows one to obtain much more information on physical 
processes with the different time and spatial intervals, which define the dynamic properties 
of the system. 

The quantum system of Bose particles serves as a physical model in theoretical descrip- 
tions both the equilibrium and non-equilibrium properties of real helium. In particular, 
many articles l|-|2l| are devoted to the hydrodynamic description of normal and superfiuid 
states of such system. A brief review of the results of the investigations within the framework 



of linear hydrodynamics has been considered in the article by Tserkovnikov 



161 . The for- 



malism of quantum molecular hydrodynamics using a mode-coupling treatment was applied 



for investigations o 
In papers 



22 



23|. 



time correlation functions, collective excitation spectrum in 
24j-|26| the theoretical approaches are proposed to the description of nonlinear 
hydrodynamic fiuctuations connected with problem of calculating the dispersion for the ki- 
netic transport coefficients and the spectrum of collective modes in the low-frequency area 
for superfluid Bose liquid. The generalized Fokker-Planck equation for the non-equilibrium 
distribution function of slow variables for quantum systems was obtained in paper by Moro- 



zov |27| . Problems of building the kinetic equation for Bose systems based on the microscopic 



approach were considered in papers 



28 



29|. For normal Bose systems, the calculations of 



the collective mode spectrum (without accounting for a thermal mode), dynamic structure 



factor, kinetic transport coefficients 



15 



16 



see the references] are carried out on the basis 



of the hydrodynamic or kinetic approaches. Nevertheless, these results are valid only in the 



hydrodynamic area (small values of wave vector k and frequency u). The papers j29|-|32| 
were devoted to the investigation of the dynamic structure factor and of collective excitations 
spectrum for superfluid helium. 

In papers HQ a generalized scheme for the theoretical description of dynamic prop- 



2 



erties of semiquantum helium has been proposed based on the method of non-equihbrium 
statistical operator. Here the set of equations of the generalized hydrodynamics is ob- 
tained and the thermal viscous model with kinetic and hydrodynamical collective modes is 
analyzed in details. The closed system of the equations for time correlation functions is 
obtained within the Markovian approximation for transport kernels. Using these equations 
the analysis of dynamic properties of semiquantum helium is carried out at two values of 
temperature above transition to a superfluid state. Similar investigations were performed in 
papers 



35 



37| for helium above a point of the phase transition. 



In general, a hard problem exists in the description of Bose systems going out from the 
hydrodynamic area to the area of intermediate values of k and a;, where the kinetic and 
hydrodynamic processes are interdependent and should be considered simultaneously. This 
is one of the urgent problems of the statistical theory of non-equilibrium transport processes 



in quantum liquid. It should be noted that in the paper by Tserkovnikov 



37| . a problem 



of building of the linearized kinetic equation for the Bose system above critical temperature 



was considered by means of the method of two-time Green functions 



40 



38 



44| in which the approach of the 



A considerable success was achieved in papers 
consistent description of kinetics and hydrodynamics of classical dense gases and fluids is 
proposed based on the method of Zubarev non-equilibrium statistical Zubarev operator 



21 



45H48|. By means of this formalism the non-equilibrium statistical operator of many- 
particles Bose system, which consistently describes kinetics and hydrodynamics, is obtained 



m papers 



49 



. l50|. The quantum non-equilibrium one-particle distribution function and the 
average value of density of interaction potential energy have been selected as parameters of 
consistent description of a non-equilibrium state. 

On the other hand, the large-scale fluctuations in a system related to the slow hydrody- 
namical processes play the essential role at the phase transition with temperature decrease. 
The construction of the kinetic equations taking into account the slow processes is the hard 
problem for the transport theory both in classical and quantum liquids. The same problem 



arises at the description of low 



tails" of correlation functions 



requency anomalies in kinetic equations related to "long 
5l|-|55|. 



The experimental investigations and model descriptions of phase transition were con- 



sidered in papers 



56 



-|60|. In this regard the studies of nonlinear second sound waves, the 



transition from laminar to turbulent flows and acoustic turbulence in superfluid helium are a 



3 



62 



71] ■ These area of science 



72 



-l76|. 



hard problems in theoretical and experimental investigations 
actively develops and is called quantum turbulence in quantum fluids 

The aim of the present paper is construction of the kinetic equations for quantum system 
taking into account the nonlinear hydrodynamic processes. The non-equilibrium statisti- 
cal operator which consistently describes both the kinetic and nonlinear hydrodynamical 
fluctuations in a quantum liquid is calculated. The coupled set of kinetic equations for 
quantum one-particle distribution function and generalized Fokker-Plank equations for the 
functional of hydrodynamical variables: densities of particle number, momentum and energy 
is obtained. Using cumulant representation a structural function of hydrodynamic fluctu- 
ations is calculated. It provides a possibility to analyzing the generalized Fokker-Planck 
equation in Gaussian and higher approximations for dynamic correlations of hydrodynamic 
variables that is important to describe the quantum turbulent processes. The generalized 
Gross-Pitaevkii equation for the order parameter which takes into account nonlinear hydro- 
dynamic fluctuations is obtained. 



II. KINETIC EQUATION FOR THE NON-EQUILIBRIUM WIGNER FUNCTION 
AND FOKKER-PLANCK EQUATION FOR DISTRIBUTION FUNCTION OF THE 
HYDRODYNAMIC VARIABLES 

Observable average values of energy density (£q)*, momentum density (-Pq)*, and particle 
numbers density (riq)* are the abbreviated description parameters at investigations of the 
hydrodynamical non-equilibrium state of the normal Bose liquid which is characterized by 
processes of the energy, momentum and masses flows. Operators for these physical quantities 
are deflned through the Klimontovich operator of the phase particle number density Tiq(p) = 




4 



where i^"^ and e^* are Fourier-components of the operators of kinetic and potential energy 
densities. Average value of the phase particles number density operator is equal to the 
non-equilibrium one-particle distribution function /i(q, p,t) = (riq(p))*, which satisfies the 
kinetic equation for quantum Bose system. 

The agreement between kinetics and hydrodynamics for dilute Bose gas does not cause 
problems because in this case the density is a small parameter. Therefore, only the quan- 
tum one-particle distribution function /i(q, p;t) can be chosen as the parameter of the 
reduced description. At transition to quantum Bose liquids, the contribution of collective 
correlations, which are described by average potential energy of interaction, is more im- 
portant than one-particle correlations connected with /i(q, p;t). From this fact it follows 
that for consistent description of kinetics and hydrodynamics of Bose liquid, the one-particle 
non-equilibrium distribution function along with the average potential energy of interaction 



49 



50|. The non- 



necessarily should be chosen as the parameters of the reduced description 
equilibrium state of such quantum system is completely described by the non-equilibrium 
statistical operator g{t) which satisfies the quantum Liouville equation: 

d 

—g{t) + iLNQit) = -e{g{t) - gg{t)) . (5) 

The infinitesimal source e in the right-hand side of this equation breaks symmetry of the 
Liouville equation with respect to t — )■ —t and selects retarded solutions (e — )■ -|-0 after 
limiting thermodynamic transition). The quasiequilibrium statistical operator gq{t) is de- 



termined from the condition of the informational entropy extremum at t 
normalization condition Sp gq{t) = 1 for fixed values (nq(p))* and {i^^^Y 



le conservation of 



49 



50|: 



gq{t) = exp I ~m - J2 ^-nitK' -EE ^-'i(P5 \ , 

I q q p J 



(6) 



where the Lagrangian multipliers /3_q(t), 7-q(p; t) are determined from the self-consistent 
conditions: 

The Massieu-Plank functional 



^t) = In Sp exp I - J] P.^{t)ef "EE ^-^(P^ ^X^P) ] 
I q q p ) 

is determined from the normalization condition. Here ((..•))* = Sp{...)g{t), ((...)) 
Sp{.. .)gq{t). 



(7) 



The system of equations ( ITT]) for the one-particle distribution function and the average 
density of potential energy is strongly nonlinear and it can be used to description both the 
strongly and weakly non-equilibrium states of the Bose system with a consistent consider- 
ation of kinetics and hydrodynamics. The description of weakly non-equilibrium processes 
was reviewed in [50]. Projecting transport equations on the values of the component of the 
vector ^(p) = ^1, p, 2m ~ 8mJ ' shall obtain the equations of nonlinear hydrodynamics, 
in which the transport processes of kinetic and potential parts of energy are described by 
two interdependent equations. Obviously, such equations of the nonlinear hydrodynamic 
processes give more opportunity to describe the process of mutual transformation of kinetic 
and potential energy in detail at investigation of non-equilibrium processes occurring in the 
system. 

In this manuscript the non-equilibrium quantum distribution function /i (q, p; t) = 
(nq(p))* is chosen as a parameter to describe of one-particle correlations. Now, to describe 
the collective processes we introduce the distribution function of hydrodynamic variables in 
a quantum system as follows: 

/(a) = / d^e-^'--\ (8) 

where a = {ctik, 0,2^, a^k} , di^ = fii^, = Pki = ik = ^k" + ^IT* Fourier- 
components of the operators of particles number, momentum and energy densities ([I])-(jl]). 
The scalar values a^k = {^ki Pk) ^k} are the corresponding collective variables. The 
operator function ([8]) is obtained in accordance with Weyl correspondence rule from the 



classical distribution function 



27| 



N 



f{a) = 6{A - a) = YlYl <^(^"^k - Omk), 

m=l k 

where A = {^ik • • • , ^TVk} are the classical dynamical variables. 

The average values /i(q, p;t) = (f;,k(p))*, f{a]t) = (/(a))* are calculated using the non- 
equilibrium statistical operator g{t), which satisfies the Liouville equation. In line with 
the idea of abbreviated description of the non-equilibrium state, the statistical operator g{t) 
must functionally depend on the quantum one-particle distribution function and distribution 
functions of the hydrodynamic variables: 



g{t) = g{... /i(q,p;t), f{a;t) . . .) 
6 



(9) 



Thus, the task is to find the solution of the Liouville equation for g{t) which has the form 
[9l) . For th at p urpose we use the method of Zubarev non-equihbrium statistical operator 



22 



45 



46 



48| . We consider the Liouville equation ([5]) with infinitely small source. The 



source correctly selects retarded solutions in accordance with the abbreviated description 
of non-equilibrium state of a system. The quasiequilibrium statistical operator gq{t) is 
determined in the usual way, from the condition of the maximum informational entropy 
functional: S{g') = —Sp{g'\iag') — Xlp 7-q(P; ^)Sp(^''^q(p)) ~ / da F{a;t)Sp{g' f{a)) with 
the normalization condition: Sp^g(t) = 1. Then the quasi-equilibrium statistical operator 
can be written as 



Qqit) = exp <^ - 5^ 5^ 7-q(p; ^)^q(p) - da F{a- 1) f{a) 

I q P 

where da — )■ {dn\f_, dP\f_, dsk}. 

The Massieu-Plank functional is determined from the normalization condition: 



(10) 



$(t) = In Sp 



exp < -^7-q(p;^)^q(p) - / daF{a;t)f{c 

I qp 



Functions 7-q(p; t) and F{a, t) are the Lagrange multipliers and can be defined from the 
self-consistency conditions: 



/i(q, P; t) = {h^{p)y = {h^{p))i, f{a- 1) = (/(a))* = {f{a))l 



(11) 



After constructing the quasiequilibrium statistical operator ( llUp . the Liouville equation 
([5]) for the operator Ag{t) = g{t) — gq(t) is written in the form: 



d 



+ iLN + e) Ag{t) 



d 



+ iLn gq{t) 



(12) 



dt ' " ' J \dt 
Time derivative of the right-hand side of this equation can be expressed through the pro- 



jection Kawasaki- Gunton operator Pq{t) 

d 



22 



27 



dt 



gg{t) = -Pq{t)iLNg{t). 



(13) 



In our case the projection operator acts arbitrary on statistical operators g' according to the 
rule 

dgq{t) 



Pq{t)g' = gq{t)^^g' + 



qp 



5(^q(p))^ 



[Sp (riq(p)^') 



[p))*Sp^'] + 



da 



dggjt) 
df{a-t) 



Sp (/(a)^') - /(a;t)Sp^' 



(14) 



Taking into account relation (fT3|) . we rewrite the equation (fT2|) as follows: 

+ (1 - PMYLn + ^ Ag{t) = -(1 - Pq{q))iLMQg{t). 
Formal solution of ( !T5|) is 



_9 

di 



where 



T,(t; t') = exp+ J - y dt' (1 - P^(t')) ilM 



t' 



is the generalized time evolution operator, that takes into account the projection. Now 
consider the action of the Liouville operator on the quasi-equilibrium operator (ITU]) : 

1 

qp 

1 





where riq(p) = iLArnq(p). Since [27 



with 



iTxa 



J{a) = (211)-^ J dxe"^^-"^ dre-'^^HL 

the second term in the right-hand side of (fTSjl can be represented: 

1 

dT{Q,{t)y{- — J{a)){Q,{t)Y-^ = 



1 

da{^F{a;t)) f dr{g,{t)y J{a){g,{t)y~^ 



da 

b 

Then the expression (fTS!) taking into account (fT8|) will take the form: 

/d 
da{—F{a;t))j{a,T)g^{t), 



where 



n^{p;r) = J dT{g,{t)yn^{p){Q,{t))-\ J{a,T) = J dT{Q,{t)y J{a){g,{t))-\ (22) 



Taking into account f l2T]) according to (fT6|) we represent the non-equihbrium statistical op- 
erator in form: 

t 

m = m + Yl j X (l-P,(f))Aq(p; r)7-q(p; 0^.^(0- (23) 

t 

j da j dt'e<''-'%{t,t'){l-Pg{t'))Jia- r)^F{a;t')g,{t') = 



t 



E/./«.-..,,™,-,„„,. 

^ — oo 

which contain non-dissipative part gq{t) and dissipative one, that consistently describe non- 
markovian kinetic and hydrodynamic processes with microscopic flows f;,q(p;r), J(a;r). 
Moreover, 

^q(p) = [nq(p), i^]- = -z^^nq(p) + (24) 

k,pi 

and we obtain the microscopic conservation law of particle density fiq. 

= [nq, H]^ = -?(q • Jq) = "^(q " Pq), (25) 

Where Jq is the flow density operator of Bose particles. Respectively, the expressions 
J(Pq;r), J(£q; r) contain the microscopic conservation law of momentum density opera- 
tor 



pet ^ 

^q = V — (p ■ q)riq(p) + -- Y\v{y)K + /^(k + q)(-A;„ + g„)]rikn_k+q (26) 

P k 



and the microscopic conservation law of total energy 



\fN ^-^ \2m 8m J m 
^ p 

VN ^ J- iy{k) - i^{-k + q) u{k) + u{-k + q) ^ 



k,a 



2V 

Non-equilibrium statistical operator is a functional of abbreviated description parameters 
/i(q, p;t) and f{a;t) contained into self-consistent conditions ( 12 ip of determination of La- 
grange multipliers 7q(p;t). The values F{a;t). /i(q, p;t) are necessary for complete de- 
scription of transport processes in system. With this aim in mind we use the ratio: 

d d ■ d / ^ ^ 

— (Tigq(p))t = ^/i(q, P; t) = (nq(p))j, ^/(a; ^) = Sp \^g{t)iLNf{a] 

Opening the average values in the right parts of these equations with non-equilibrium sta- 
tistical operator (123|) we obtain the system of transport equations: 



o^/i(q, p; t) + i ^ /i(q, p; ^) = Yl ^^^^pi,p-^ ~ '^pi,p+|)^2(p, q, pi, k; t) + 
^ k'pi' 



k'pi'j 



J2 I ^it'e^^*'-*)y,„„(q, q', p, p'; t, t')7-q'(p'; ^0 " (28) 

) 

t 

J da J c/tV(*'-*Vnj(q,P,;t,t')^i^(a; t'), 



p{a;t) + ^{J{a))l = "E / ^^V^*'"*) |-^,„(a, q', p'; t, t')7-q'(p'; ^0 + (29) 

t 

I da' I dt'e<''''^^^jj{a,a'-t,t')-^F{a'-t'), 

— oo 

where 

5'2(q,P, Pi,k;t) = Sp(a+ k-q^q(p)a k^^'/(i)) (30) 
is the pair quasi-equilibrium distribution function of Bose-particles, i/(k) = ^^z/(k), 
'/'nn(q) P) q') p'; ^) i^') is the transport kernel (memory function), which describes dissipation 
of kinetic processes: 

¥'nn(q,q',P,p';t,0 = {in{ci,p-t)f,{t,t')Uc^,p'-t',T)Y^. (31) 

10 



Respectively (/?„j(q, p, a'; t, t'), Lpjn{a, q, p; t, t') are the matrix with elements 

(/;„j,(q,p,a';t,0 = (4(q, p; t)T;(t, t', t));', (32) 

(/?j,„(a,q,p;t,t') = (Jj^a; t)Tg(t, t')4(q, P; t', t"))*'. 

These elements are the transport kernels, which describe dissipation between kinetic and 
hydrodynamic processes. Lpjj{a, a'; t, t') is the matrix with elements 

ipj,j^{a,a']t,t') = {ij^{a;t)fg{t,t')ij^{a;t',T)Yg, (33) 

which describe dissipation of hydrodynamical processes in quantum Bose fluid. The trans- 
port kernels (pT]) - fl33p are constructed on generalized flows: 

/„(q,p;t)= (l-F(t)) Aq(p), (34) 

/j,(a;t) = (l-P(t)) JKa), (35) 

where P{t) is the generalized projection Mori operator, which is constructed on the operators 
hq{p), f{a) and corresponds to structure of projection Kawasaki- Gunton operator Pq{t) 
(fl^ . It is signiflcant that transport kernel contain contribution from quantum diffusion in 
coordinate and momentum space and contribution from generalized function "force-force". 
One can readily see by substituting (IMl) into (^„„(q, q', p, p'; t, t') and open contribution 
from kinetic and potential parts of ( IMl) : 

fnn{q, p, q, p; t, t') = ^q ■ L>„„(q, p, q, p; t, f) ■ q - — q ■ D„F(q, p, q , p; t, t') - 

DFn{(i, p, q', p'; t, t')— q' + ^^^(q, p, q', p'; t, t'), (36) 

m 

where 

Dnniq, P, q', p'; t, t') = ((1 - P{t)) ■ pnq(p) • f,(t, t') (l - P{t')) ■ p'n^{p% (37) 

is the generalized diffusion coefficient of quantum particles in the space {q, p}. Other 
components have the following structure: 

DnFici, P, q', p'; t, t') = ((1 - Pit)) ■ pnq(p) • f,(t, t') (1 - Pit')) ■ F^,ip'))l, 
DpnicL, P, q', P'; t, t') = ( (1 - Pit)) ■ Fq(p) ■ f^it, t') (1 - Pit')) ■ p'nq.(p'));, 
DM<i, P, q', P'; t, t') = ((1 - Pit)) ■ Fq(p) ■ f,(t, t') (1 - Pit')) ■ (p'));, 

11 



where 

'N 



IV-f^ y / J \ )\ Pi,p-2 Pi,P+2^ Pi,p-f IVf^ p^^p 

k,Pi 



For a detailed study of the mutual influence of kinetic and hydrodynamic processes we 
will allocate the "kinetic" part in quasi-equilibrium statistical operator using the operator 
representation: 

1 

m = gf"{t) - I daF{a; t) J drQiF\r){g';-{t)y f{a){g'^^^ (38) 



where 

g'^^it) = exp{ - - 5^7_q(p;t)nq(p)} (39) 

q,p 

is the quasi-equilibrium statistical operator, which is the basis of kinetic level of description. 
The value Q{F\t) satisfies the operator equation: 

1 

Q{F\t) = 1- J daF{a;t) J dTQ{F\T){gf^{t)r f{a){g^^-{t))-\ (40) 



One use the representation (138|) for determining of Lagrange multiplier F{a\ t) from the 
self-consistent condition: 



/(a; t) = {f{a))l = {f{a))l^ - j da'W{a, a'; t, T)F{a'- 1), (41) 

where 

1 

W{a,a'-t,T) = I dT{f{a)Q{F\r){g':;-{t)Yf{<^){^^^^ (42) 



is the structural function, in which the averaging implement with quasi-equilibrium statis- 
tical operator f HOj) . From (HT]) we find F{a\t) in form: 

F(a; = - y da's f {a'; t)W~\a, a'; t), (43) 

Function W~^{a, a'; t) is the solution of integral equation: 

j da"W{a, a"; t)W.i{a" , a'; t) = S{a - a'), (44) 

and 

Sf{a- 1) = f{a- 1) - (/(a))L = (/(a))* - (/(a))L (45) 



12 



are the fluctuations of distribution functions of hydrodynamic variables determining as dif- 
ference between the complete distribution function and one averaged with operator Qq^"'{t). 
Taking into account ( H3|) the quasi-equilibrium statistical operator Qgit) can be presented 

as: 

1 

m = Q'rit) + Jd(^J da"W.,ia\a;t) J rfrg(F|T)(^f "(t)) '^(^))'-^ (46) 



or as: 

1 

= "(t) + jdaj da" j dTQ{F\T)W.i{a',a; t)f{a;r)6f{a';t)g^'''{t), 



where 

f{a;T) = {g';-{t)Yf{a) {gf-{t))-\ 

An important point is that expression (H3|) is equation for determining F{a;t), because 
the function W^i{a, a']t) depends on F{a;t) according to structural function W{a,a';t,T), 
which in turn depends on Q{F\t) ( HOl) . Taking into account (H3l) and ( H71) we rewrite the 
equation system (128|) . (129|) in form: 



d 



/i(q,p;t) + 2-^^^— ^/i(q,p;t) = ^ f>k((^p,,p_k - 5p,,p+i£)c/2(p, q, Pi, k; t) + 



9t ' ' m 

k'pi' 



t 



J2 j rft'e^^*'-*Vn„(q, q, P, P; t, t')7-q'(P ; + (47) 

P'-oo 

t 

j da' j da" j rftV(*'-*Vnj(q,P,;t,t')^W^-i(«',a";0'5/(«;0, 



— (5/(a; t) + — / (ia' / da"v(a, a";t)W^i(a" , a;t)6f(a';t) 
at da J J 

t 

^ j dtV(*'-*)^¥,,„(a, q', p'; t, t07-q'(P ; ^0 - 



(48) 



q p 

t 



da' I da" j dt'e<''-'^^^jj{a,a';t,t')-^W.,{a',a";t')6fia";t'), 



where was taken into consideration ■§^{J{o))\in — ^ §t^f {'^)Ykin-> 
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via,a';t)= I da" I drSp ( J(a)g(F|T)/(a"; r)^^'*"(t) ) W^_i(a", a'; t) 



(49) 



The two limiting cases follow from the system of transport equations HTj), (HSj) . First, 
unless we consider nonlinear hydrodynamic correlations then obtain a kinetic equation for 
Wigner function of quantum Bose particles. Second, if we take no account the kinetic pro- 



cesses then obtain a Fokker-P 
to results of Morozov article 



ank equation for distribution function /(a; t), that corresponds 



^/(a; ^) + ^y" dav{a,a)f{a;t) = 

t 

j j da'K{a,o!-t-t')^ j da"W_i{a',a")f{a";t), 



(50) 



where fH9|) transforms to v{a,a') from |25|: 

v{a,a') = j da"Sp J(a)/(a")W^-i(a'.a") 
The function 

W-i{a" , a) = W-i{a)[S{a" — a) + r{a" , a)] 
is the solution of integral equation 

da"W{a, a")W{a" , a) = 6{a — a) 

and contains the singular and regular components, 

W{a,a') = Sp (f{a)f{a' 

is the structural function 



(51) 



(52) 



(53) 



(54) 



25|. Accordingly K{a, a'; t) is the matrix with elements: 
Kif{a, a'- 1) = Sp (Xi{a)f,{t, t')Xf{a')) , (55) 



where Xi{a) = (1 — P)Ji{a) are the generalized flows with projection operator PA 

i(a, a')Sp(A/(a')). If we neglect in system of transport equation 
(149|) by memory effects on hydrodynamical level then obtain the equation system: 



d 



^/i(q, P; t) + i-^^^fi{q, p;t) = ^k(5p,,p_k - p^k)^2(p, q, Pi, k; t) + (56) 



dt 

t 



kpi 



^ j dt'e'^^' *Vnn(q,q',P,pV,^')7-q'(p';^') - y c/aVnj(q,P,;a')^W/_i(a';t)5/(a;t), 
q' p'-oo 
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^5/(a; t) + -^via; t)Sf{a; t) = - ^ -^ipjnia, q', p')7-q'(p'; t') - 

q' p' 
d d 

in which the memory effects on kinetic level are preserved. In this system the following 
designations are used: 



(/7„j,(q,p,a) = j dte^*(t(q,p)f;(t)/jXa))L, (57) 

— oo 



^JiJf{a) = I da' I dte''{ij,{a)f,{t)ij^{a'))l^. (58) 

— oo 

v{a; t) is the contribution from singular part of generalized velocity v{a, a'; t) = v{a; t)6{a — 
a') + u{a, a', t), during which v{a; t) = {J{a, t))\^^. An additional point to emphasize is that 
only contribution from singular part of structure function W-i{a,d\t) is in this equation 
system, namely W^i{a\t). Such local approximation can be applied near a critical point 
when the values, that strongly fluctuate, are the hydrodynamical variables and the long- 
wave components of order parameter. In particular, for superfiuid helium it can be the 
averaged values of wave function (^q)o, that describes the Bose-condensation phenomenon. 
These value can be taken into account if to generalize the quasi-statistical operator ( !T0|) 

as 



Qqit) = exp \ -$(t) - 5^ 5^ 7-q(p; ^)^q(p) " &q(i)^q- / dtt F {a; t) f{a) I 

I q P q ^ 



(59) 



where the parameter 6q(t) is determined from self-consistent condition: 

(^q)' = (^q)^ (60) 

Consistency with the structure of quasi-equilibrium statistical operator gq{t) the kinetic 
equation for (^q)* is presented in form 



t 



|(^q)* = (^q)* + 5^ I dtV(*'-*)y,^^(q,q';t,t')&-q'(0+ (61) 

l' -oo 

t 

J2 I rftV(*'-*VM'n(q,q',p;t,t')7-q'(p';0 + 

q'P'-oo 

t 

j da' j do!' j dt'e^'-''-''>ip^j{fi;t,t')-^W_i{a',a";t')5f{a;t'), 
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where = zI/Ar^q, (f^q,{c[, q'; t, t'), (y9^„(q, q', p'; t, t'), if^j{q;, t, t') are the transport kernels 
which describe the dissipative processes in a system. 

The equation system (H7|) and (H5|) is modified also depending on flS^ . The equation 
(j6T|) can be called as Gross-Pitaevskii equation \r\ . \l\ which takes into account the dissi- 
pative kinetic and hydrodynamic processes. This formulation deserves of separate detailed 
consideration. 

Discussion 

The non-equilibrium statistical operator of consistent description of kinetic and non- 
linear hydrodynamical fiuctuation for quantum Bose system is obtained. For consideration 
of kinetic processes the non-equilibrium one-particle Wigner function is used as parameter 
for abbreviated description. The distribution function of hydrodynamical variables is chosen 
for study of non-linear hydrodynamic fiuctuations. The equation of type Fokker-Plank one, 
which coupled with kinetic equation, is obtained. The non-equilibrium distribution function 
provides a possibility to calculate the averaged values {{cii . . . ai)Y = j da {ai . . . di) f{a; t) 
and particularly, for example, to obtain for fiux density operators Jq a chain of type Reynolds 
equations for values (Jq)*, (Jq Jq')*, (JqJq'Jq")* as in classical case [79]. It is important to 
research of quantum turbulence phenomena. 

A hard problem at the examination of nonlinear fiuctuations based on equation system 
(H7|) . f l48p is the calculations of structural function W{a; t) and generalized hydrodynamical 
velocities v{a]t). One consider the function W{a;t) and rewrite as 

1 

W{a-t) = (/(a,r))L = ^/ dxe'^^^ J drSp((^J-(t))^e-^(^J-(t))i--). (62) 



Now using the cumulant expansion this function represent in form: 

W{a; t) = - — - I dxe " "-"i , (63) 



(27r)5 , 

where Yl, ^aO-a = a^l^q + X2 ■ Pq + X^E^, 



1 



M^{t) = (a„(r))L = J rfrSp((^J-(t))-a„(^J-(t))^--), (64) 



M^aAr) = {da{T)do,,{T))ii^ - (a„(r))*fci„(a«, (t))^^^ (65) 
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are the cumulant coefficients averaged with quasi-equihbrium statistical operator ^^'"(t). 
Such calculation of structural function (163|) enables one to consider the Fokker-Plank equa- 
tion in Gauss approximations and higher in cumulant averages and by this means to obtain a 
chain of type Reynolds equations for time correlation functions ((a/...aj))* in correspondence 
approximations. It is important to description of nonlinear fluctuations. The generalized 
lydrodynamical velocities v{a\ t) will be calculated in the same manner as for classical case 



79 



83 



84| through the cumulant representations in Gauss approximations or higher. Then 



we can present the transport kernels in the Fokker-Plank equation in type mode-coupling 



form 



85| as in classical case. For quantum statistics its invites especial investigation that 



will be done in next papers. 
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